Nuclear vorticity and the low-energy nuclear response - Towards 

the neutron drip line ^ 



> 
O 



p. Papakonstantinou"''', J. Wambach", E. Mavrommatis'', V.Yu. Ponomarcv"'^ 

Institut fiir Kernphysik, Technische Universitdt Darmstadt, Schlossgartenstr.9, 
D-64289 Darmstadt, Germany 
^Physics Department, Nuclear and Particle Physics Section, University of Athens, 

GR-15771 Athens, Greece 

Abstract 

^«~^ [ The transition density and current provide valuable insight into the nature of 

■ nuclear vibrations. Nuclear vorticity is a quantity related to the transverse tran- 

' sition current. In this work, we study the evolution of the strength distribution, 

related to density fluctuations, and the vorticity strength distribution, as the neu- 
tron drip line is approached. Our results on the isoscalar, natural-parity multipole 
response of Ni isotopes, obtained by using a self-consistent Skyrme-Hartree-Fock 
-|- Continuum RPA model, indicate that, close to the drip line, the low-energy 
\^ ' response is dominated by L > 1 vortical transitions. 

Keywords: Zero sound, nuclear vorticity, exotic nuclei, collective excitations, self- 
consistent continuum-RPA. 
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2 '. 1 Introduction 

, One of the aims of the ongoing theoretical and experimental studies devoted to nuclei far 

' from stability is to illuminate the effects of the excess protons or neutrons on the nuclear 

structure and response. On the neutron-rich side of the nuclear chart, in particular, a 
wealth of exotic phenomena is anticipated, due to the formation of neutron skins and 
halos, the large difference between the Fermi energies of protons and neutrons and a 
possible change of shell structure. 
^ ] Electron scattering experiments with radioactive beams are currently being dis- 

cussed, opening the possibility of detecting, besides the transition charge distributions, 
also transition current distributions in a broad region of the nuclear chart. The tran- 
' sition current density (TC), studied in conjunction with the transition density (TD), 

M I provides valuable insight into the nature of nuclear vibrations, in particular their "zero- 

sound" nature. Theoretically, the TD and TC associated with a particular type of 
electric excitation, approach ideal hydrodynamical behaviour in the case of very collec- 
tive states such as giant resonances (GRs). It is predicted P that the properties of 
GRs in drip-line nuclei are close to the ones of /3— stable nuclei in this sense. At lower 
energies, though, the behaviour of the TD and TC may deviate significantly from this 
picture, as several examples demonstrate. Microscopic studies of convection currents in 
stable nuclei predict that a variety of flow patterns is possible [2113 El; for instance, 
the isoscalar quadrupole GR of ^''^Pb is nearly irrotational, whereas the correspond- 
ing low-lying state exhibits considerable vorticity. In (stable and unstable) neutron-rich 
nuclei, "pygmy" (SJ |B] or "soft" [3 HI |5I dipole modes are observed, interpreted as os- 
cillations of the excess neutrons against the isospin-saturated core; also, toroidal dipole 
modes are predicted to exist at a somewhat higher energy E3] . A special feature 
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of the response of very neutron-rich nuclei is the so-called threshold strength, i.e. the 
considerable amount of strength predicted to occur just above the neutron threshold; 
it is of isoscalar and non-collective nature and characterized by spatially extended TDs 
and TCs, owing to the loosely bound neutrons El El El E]- 

The structure of the particle continuum and the nature of low-lying strength in 
unstable nuclei are of particular interest, as they provide information on decay properties, 
polarizabilities, shell structure etc. Close to the drip line, collective "low-lying" states 
such as 2+ may be located in the continuum. The energy of 2+ states as well as other 
Ohio modes, may be very low even for nuclei described traditionally as closed-subshell 
configurations, if the shell structure tends to "melt" at the extremes of isospin Such 
transitions can hardly be expected to be irrotational. 

A standard procedure in studies of nuclear convection currents is to define a velocity 
field, as dictated by hydrodynamical considerations. In Ref. it was shown that the 
degree to which a transition is irrotational, can be quantified in terms of the nuclear 
vorticity density. The latter is directly related to the transverse component of the TC and 
therefore, together with the TD, it completely specifies a transition and unambiguously 
reflects the zero-sound character of nuclear vibrations. 

In this work, we examine the response of the isotopes 56,78,ii0]\jj ^ ranging from 
N — Z and close-to-stable {N = 28) to extremely neutron-rich {N — 82). We focus 
on the flow properties of the various transitions. The isoscalar strength distributions 
Sp{E), related to density fluctuations, and respective vorticity strength distributions 
Sui{E) of multipolarity up to i = 4 are calculated, using a self-consistent Skyrme- 
Hartree-Fock plus Continuum-RPA method. We consider the vorticity related to the 
convection current. The evolution of the distributions as the neutron number increases 
is discussed. 



2 Definitions and method of calculation 

Let us consider the response of spherical, closed- (sub) shell nuclei to an isoscalar (IS) 
external field VLAiir) = r^YLM{0, </>) {L > 0). The IS TD Sp{r) and TC J(f) character- 
izing the excited natural-parity state \LM) of energy are determined by their radial 
components, SpL{r) and jLL±i{r) respectively |21E|. From j{r) and the ground-state 
density p{r), a velocity fleld can be deflned, {[(r) = j{r)/p{r). Under the assumption 
that a single transition exhausts the available strength of an external potential V{f), it 
can be shown that ti(r) oc VF(7^ ^Hj) which implies irrotational flow, V x u = 0. In 
classical hydrodynamics, the flow is always irrotational in ideal, non-viscous fluids. 

The IS TD and TC are related through the continuity equation (CE) '2', which we 
write here for the radial components: 



E, L d L-l.. ,. L + l^d L + 2.. 

yJp"^^'^ = V ^^^d^ " ^W-i(0 - V^^(d^ + ^ WiM- (1) 

The transverse part of the current, j'', with radial components jLL±ii''')^ does not 
contribute to the CE since, by definition, V • j'*' = 0. This equation does not define 
j''' uniquely; following Ref. TB", one may set jj^L+ii^) ~ jLL+iir)- Then ji;i^_i{r) is 
determined from V-j*"^ = 0. The remaining part of the current is related to SpL through 
the CE. Therefore, measurements of Sp^ and jLi_|_i(r), as demonstrated for example in 
the analyses of inelastic-electron-scattering data in Ref. determine the transition 
completely. 

The vorticity density LJ(r) associated with a transition from the spherically symmetric 
ground state to the excited state \LA'I), as defined in Ref. jj, equals the curl of the 
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transverse current introduced above |17| . 

uj{r) = V X j'^ir) EE (2L + l)-'/'cuLL{r)Yli{f). (2) 

The radial part LULL{r) is given by 

/2L + 1, d L + 2, 
coLLir) = y — + ^)jLL+i(r). (3) 

The r^— moment of ujll is identicaUy zero; its r^+^— moment may be used as a measure 
of its strength |4j. For irrotational and incompressible flow, e.g. in the Tassie model, 
LULL vanishes. 

The strength distribution Sp{E) = J2f \{f\VLM\0)\^S{E - Ef), where |/) are the 
final states excited by the field Vlm and Ef their excitation energies, is related to 
the r^— moment of 5pl{t). Since we are dealing with continuous distributions, we 
write the strength in a small energy interval of width A£^ at energy E as Sp{E) = 
[/q drSpL{r)r^~^'^]'^ / AE. We define the vorticity strength distribution S^{E) in a sim- 
ilar way, namely S^{E) — drcjLi(r)r^+'']^/A_E. Comparison of Sp{E) and St^{E) 
allows a discrimination between collective states of longitudinal zero-sound character 
and excitation modes which are highly vortical and possibly of transverse zero-sound 
character. The latter case applies, for example, to the toroidal dipole mode I10| . 
Note that, in order to quantify how "collective" (in terms of S'^^) a seemingly strong 
transition is, and to what extent it is a transverse-zero-sound candidate, one would need 
to establish appropriate sum rules for S^j. 

The quantities introduced above are calculated using a Skyrme - Hartree-Fock (HF) 
plus Continuum - RPA model. For the HF ground - state, the numerical code of P.- 
G. Reinhard ^H] is used. The calculation of the response function is formulated in 
coordinate space, as described in |19ll2()ll^l22| . First, for given energy E, multipolarity 
L and isospin charachter r^, the radial part G^^^^ (r, r'; E) of the unperturbed ph Green 
function G^{E) is calculated. The indices pi — 1,2,... enumerate the operators whose 
propagation is described by the Green function, namely 

Ylm, [Yl®{V^ +V'^)]lm, [Fiii (V ± V')]la/. 

Spin-dependent terms have been omitted in the present calculation. The continuum is 
fully taken into account, as described in |19[ I21| . A small but finite Im£^ ensures that 
bound transitions acquire a finite width T^. The RPA ph Green function is obtained 
by solving the equation 

G«PA(ii;) = [1 + G\E)V,,,]-^G\E) (4) 

as a matrix equation in coordinate space (represented by a radial mesh) , isospin character 
and operators pn. The ph residual interaction Vrcs is derived self-consistently from the 
Skyrme-HF energy functional The radial functions (5pL, Jll-i and Jll+i, at 

given energy E, can be calculated from the RPA Green function; then, loll is obtained 
using Eq. Q ElESl. 

3 Results and discussion 

We have performed calculations for the IS natural-parity response of Ni isotopes for 
L < 4. Since pairing is not included in our calculation scheme, only the doubly closed 
isotopes ^''Ni {N = 28), ^*Ni (TV = 50) and "°Ni {N = 82) have been studied. 
The nucleus ^^Ni lies close to the valley of stability, while ^^°Ni serves as a numerical 
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example of a nucleus with extreme neutron excess lying close to the drip line. In all 
cases we examined, we have checked that the CE, Eq. is fulfilled to a satisfactory 
degree. We present results for the IS 2+, 3~ and 4+ strength distribution Sp{E) and the 
corresponding vorticity strength distribution S,^{E). The distributions are normalized 
to 1, i.e. Sp^uiiE)/ Sp^c^^ tot is shown, where Sp^uj, tot is the integrated strength in the 
region 0-50 MeV. For the results shown we have used the Skyrme-force parametrization 
SkM* 26 , tailored to describe GRs of stable nuclei and employed in previous studies of 
the response of exotic nuclei as well (eg. in [TI ITTIIT^IOI ). 

As we observe in Fig. 1, most of the IS quadrupole strength is concentrated in two 
peaks, a low-energy collective state and the GR. The width of the GR increases and its 
energy decreases as we approach the neutron drip line. The GR carries a large fraction 
of S'p, tot and little vorticity. It has been verified that, in the three Ni isotopes studied, 
its velocity field resembles that of ideal hydrodynamical flow. On the contrary, the 
lowest 2+ state is strongly vortical. The octupole response of ^^Ni (^*Ni ), Fig. 2, is 
dominated by one (two) collective state(s) carrying almost no vorticity strength; a large 
part of the Sui lies in the 3hoj region. The situation is quite different in the case of ^^"Ni 
, where highly vortical threshold strength appears at low energy. In order to illustrate 
the difference, in Fig. 3 we plot the velocity field corresponding to the 9.3 MeV state 
in ^^Ni and the 3.4 MeV transition in ^^°Ni , in an arbitrary scale. For ^^^'^^Ni , it 
is interesting to notice a very weak octupole state with considerable vorticity strength, 
at around 15 MeV. The hexadecapole strength distribution of ^^-'^^Ni , Fig. 4, appears 
fragmented. Vorticity strength is carried by all peaks. In ^^°Ni a striking amount of 
vorticity develops at very low energy. 

The neutron-rich nucleus ^^Ni is possibly doubly magic Q7^. We notice that its low- 
energy response does not differ dramatically from the response of ^^Ni . The transition 
strength of ^^'^Ni is located in the continuum in all examined cases. We have performed 
calculations also for the 0^ and 1~ response of the three isotopes. In ^^°Ni , the threshold 
energy for these multipoles is higher than for the L = 2,3,4 multipoles. The threshold 
strength of ^^°Ni for L > carries increased vortical strength. This result should be 
anticipated given the single - particle character of the threshold strength. The present 
results indicate that, close to the neutron drip line, the low-energy response is dominated 
by L > 1 transitions of transverse character. 1~ as well as 3~ excitations may gain 
importance with respect to 2+ and 4+ in the special case of nuclei described as n^— closed 
configurations. 

We have also used the Skyrme force MSk7 [2HI, whose parameters were determined 
by fitting the values of nuclear masses, calculated using the HF+BCS method, to the 
measured ones, for 1888 nuclei with various values of |iV — Z\/A. The two forces have 
similar nuclear-matter properties, except for the effective mass m* . The SkM* force 
{m* /m = 0.786) predicts that the last occupied neutron state of ^^°Ni is bound. The 
nucleus appears soft against excitations of vortical nature. In the case of the MSkT 
force [m* /m = 1.05), the last neutron state of ^^"Ni appears as a resonant state at 
positive energy and therefore the nucleus lies beyond the neutron drip line. In general, 
the results obtained with the two Skyrme parametrizations agree qualitatively. The 
strength distributions are systematically shifted towards lower energies when MSk7 is 
used. 

We should stress that the nucleus ^^''Ni was studied as an academic example of 
a closed-shell Ni isotope close to the neutron drip line, as has been done before in 
Ref. Illj . This particular choice was dictated by the limitations of our model, which 
does not account for pairing correlations. According to SHF results, the iV = 82 closure 
may still be valid in the Ni region TT, although not conclusively. Ideally, one should 
perform a (preferably self-consistent) Quasiparticle-RPA (QRPA) calculation including 
the full continuum, as proposed in Ref. 29 . The correct treatment of the continuum 
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when pairing is present, in the ground state as wcU as in the excited states, is particularly 
important in weakly bound neutron systems, where the pairing field couples the neutrons 
in the bound states with those in the low-energy continuum [301121 EH- In general, 
pairing correlations are not expected to modify dramatically the electric transitions of 
even-even nuclei (the same may not hold for magnetic transitions, however) . They 
can introduce new transitions and change the energy and strength of the ones that we 
have already examined |^|^. Therefore, we expect that our conclusions would still 
hold qualitatively, after pairing is taken into account. Of course, they are restricted to 
spherical nuclei. 

In conclusion, our results for the IS, natural-parity response of Ni isotopes, obtained 
using a self-consistent Skyrme-Hartree-Fock -|- Continuum RPA model, indicate that, 
close to the neutron drip line, instabilities develop in the form of L > 1 transitions 
of transverse character. We have used the vorticity strength distribution to quantify 
the irrotational character of the various transitions. Future work should account for the 
spin current, besides the convection current. It would be useful to examine the magnetic 
response as well, for a more complete picture. Ideally, a chain of even— iV isotopes should 
be studied using a QRPA method including the continuum. The proton drip line is also 
interesting to explore; modifications of the strength distributions, in the direction of 
sharper low-lying transitions, are expected due to the Coulomb barrier. 
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Figure 1: Quadrupole strength distribution (left) and vorticity-strength distribution 
(right), norniahzcd to 1, for the isotopes ^^''^^'i^^Ni . The Skyrme parametrization 
SkM* has been used. 
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Figure 2: As in Fig. 1, octupole strength distributions. 
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Figure 3: Left (Right): Velocity field for the octupole excitation of "''^Ni at E = 9.3 MeV 
("°Ni at E ^ 3.4 MeV) plotted up to a radial distance of 6 fm (7.6 fm). The scale of 
the velocity amplitude is arbitrary. The axis is along the page. 
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Figure 4: As in Fig. 1, hexadecapole strength distributions. 
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